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Abstract 

The accuracy of the vacuum saturation hypothesis is discussed using the examples of 
vacuum expectation values of four-quark operators and the parameter B, which deter- 
mines the short-distance contribution to the K° — K° mixing. 
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1. - The standard model of strong and electroweak interactions predicts the existence of 
transitions changing the strangeness by two units. These processes are due to exchange of 
intermediate vector bosons and lead to a possibility of K° — K° mixing. A calculation of the 
corresponding effective Hamiltonian on the basis of the standard model was first performed 
by Gaillard and Lee fl[] by expanding the interaction in the inverse mass of the intermediate 
boson. The comparison of the result of this calculation with the experimentally observed mass 
difference of Kl and K$ mesons led to a prediction for the c-quark mass before experimental 
detection of J/ip particle. In spite of the considerable theoretical uncertainty, involved in 
evaluation of the matrix element (K°\H e ff\K°) the predicted quark mass happened to be in 
rather good agreement with the measured one. However, the value of the Kl — K$ mass 
difference is quite sensitive to the numerical values of such parameters of the standard model 
as the masses of the t- and c-quarks, the mixing angles between various quark generations and 
the general structure of the model. Since the pioneer calculation by Gaillard and Lee, a good 
deal of work has been carried out to refine its results. 

Account was taken of the strong-interaction corrections to H e ff in the leading logs approx- 
imation and for the three-quark generations by Gillman and Wise || . The obtained expression 
for the Kl — K$ mass difference is of the form 



Am = m L — m s = 2Re M 12 



where 



^^l A ^^ + Ai ^ + 2AcA ^ ln ^H^P)J 2m K 



Here the standard notations were used: Aj = VidV*., i = u,c,t being the elements of the 
Kobayashi-Maskawa matrix; rj^s are some coefficients due to strong- interaction contributions, 
O is the four-quark operator O = (sl7^l) 2 with AS = 2. 

Equation (1) was obtained with the help of the Wilson expansion at short distances, so 
this contribution to the mass difference Am should be more correctly called the contribution 
of short distances. Actually, the existence of the other contributions to the mass difference was 
demonstrated by Wolfenstein ||. These new contributions are essentially different from those 
under discussion and come from long distances. We shall not dwell upon them and reserve the 
term "mass difference" and notation Am for the short- distance contribution only. 

The strong interactions enter eq. (1) through Wilson's coefficient rj^s, which can be reliably 
calculated because of the property of the asymptotic freedom of QCD, and through the matrix 
element of the four-quark operator O, being responsible for the AS = 2 transitions. It is this 
matrix element we are going to discuss. It should be stressed that this problem is of gravely 
nonperturbative nature and, thus, it is a hard nut to crack by any of the existing techniques. 
However, because of the great importance of the numerical value of the matrix element several 
attempts have been undertaken to evaluate it. 

Gaillard and Lee have estimated this matrix element by means of the vacuum saturation 
hypothesis with the result 
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fx being the K- meson decay constant. It is convenient to parameterized the exact matrix 
element with the help of dimensionless parameter B as follows: 

(K°\d\K°) = B(K°\d\K°) vs . (2) 

A lot of efforts has been made to bypass the vacuum saturation hypothesis and to calculate 
the parameter B independently, the obtained results diverge from each other [4-12]. 

This work aims at giving an analysis of attempts to compute B within the QCD sum rules 
approach. 

First, Chetyrkin et al. |l2j applied the technique of the QCD sum rules to the three-point 
correlator comprising the operator O and the interpolating currents of the i^-mesons with the 
result B = 1.2 ± 0.1. Second, after completion of that work there appeared similar in its spirit 
the calculation by Pich and Rafael [10]. They used the two-point correlator of the four-quark 
operators and obtained B = 0.38 ± 0.09. 

It is of interest to trace the cause of this marked difference since the sum rules method, as 
a rule, works with an accuracy of about (20 -v- 30)%. 

The work is organized as follows. In sect. 2 we briefly describe the calculation of B by 
means of a three-point correlator and discuss the accuracy of the vacuum saturation hypothesis 
employing there to estimate the vacuum expectation values (v.e.v.) of four-quark operators. 
Section 3 is devoted to a comparison of this calculation with the one performed by Pich and 
Rafael. Section 4 contains our conclusions. 



2. - The starting point of the work [12| is the following representation: 



(K \O\K°) = lim lim (p 2 -m 2 K )[K 2 -m 2 K )p^K u (3) 



Tn.KfK? / dxdy(0\TfJx)jl(y)O(0)\Q)exp\ipx-iKy], 



which can be obtained with the reduction formulae. Thus, to find the matrix element (3) one 
needs to compute the function 

Tp,(p, q) = i 2 J dxdy exp[ipx - iqy){Tj*(x)d(y)jl(0)} = (4) 

= p^q v T{p 2 , (p — q) 2 , q 2 ) + other structures, 
where = dj^j^s is the interpolating field of the i^°-meson: 

(o\ 3 l(o)\K (p)) = tpj K , f K = i.nu 

The function (4) can be reliably calculated at small q and large Euclidean p 2 > 1 GeV 2 by 
means of the (somewhat modified) operator product expansion technique fI3| . On the other 
hand, the matrix element (2) is connected with the amplitude T(— t, — 1,0) = T(t) by the 
dispersion relation 

r 7 p( s ) i • , 2 (^°|o|^ ) a 

T(t) = / ds^-^- - subtractions = fl\ — 1 ' ' + =- + . . . , 5 

J s + t it + m 2 K ) 2 t + m 2 K 
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where the one-pole contribution corresponds to transitions of the i^°-meson to other (different 
from K°) states and dots stand for the higher-state contributions. 

Within the vacuum saturation approximation the function T liv {p,q) assumes the form 

8 /" 
T lv = 3 n Ma (p)n va (p - q), n^(p) = i J dx exp[ipx}(Tjl(x)s L (0)-f a d L (0)} (6) 

and the resulting value of B proves to be B vs = 1. Thus, there remains to compute only 
the function A^ u = — TYf, which is responsible for all the departures from the vacuum 
saturation prediction for B. 

The computed result for A^ v is (only local operators with dimension < 6 were taken into 
account) |12|| 

\u(p, q) = Ptfv |- 5 (^)^^F " Mdssd) - A(ddss} + (7) 

+2(ssss) + 2(dddd) H s{gdG tlu a tlu d) 1 other structures, 

v ' N ' 247r 2 J p 2 (p - q) 2 

where the designation (dssd) stands for {dL"i a SLSLl a dL) and so on; G 2 = G® U G® U ; G^ v = G® u t a ; 
tr(t a t b ) = {l/2)8 ab ; = W)h*M- 



On putting q = and using the method of finite energy sum rules fllj] to connect the 
phenomenological (eq. (5)) and the theoretical (eq. (7)) representations of the function T(t), 
one gets 

-f K m\{B - 1) = / p th (s)(s + m 2 K )ds = (247r 2 )- 1 m s (gda^G^d) 

o JO 

or 

m s {gda flv G tlv d) 



B-l 



\<QK 2 f^m 2 K 



To estimate the v.e.v.'s of four-quark operators appearing in eq. (7) the hypothesis of vacuum 
saturation has been used, which leads to the vanishing of each of these v.e.v.'s. Assuming the 



relation [15 

(gda^G^d) = m 2 {dd) , m 2 = (0.8 ± 0.4) GeV 2 , 
one finds for the renormalization group invariant quantity 

5 = 5(/i)(a s (/i))- 2/9 = 1.2±0.1. (9) 



This estimate is in good agreement with the hypothesis of vacuum saturation; within this 
approach this means the absence of operators giving appreciable contributions to A^. 

2 Having in mind a solid theoretical basis of our calculation we have just updated the prediction for Bk by 
using present values of relevant parameters. The quantity B — 1 in eq. (8) is expressed through well known 
parameters that did not change much during last years. The bulk of the change of the parameter B is due to 
the normalization (the change in the factor a s (1.2 GeV 2 ) -2 / 9 in eq. (9)). The updated version of eq. (9) with 
a s (1.2 GeV 2 ) = 0.69 reads now: 

B = 1.0±0.1 
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In principle there may exist three kinds of extra corrections to A^: 

i) corrections of higher order in a s ; 

ii) corrections due to bilocal operator; 

iii) corrections due to local operators. 

Let us consider them in turn. The first contribution is suppressed by the factor a s /7i ~ 0.1 
and can hardly change the result (9) to a large extent. The second one is suppressed by an extra 
power of p~ 2 and does not contribute at all within the finite energy sum rules approach. (This 
statement is valid as long as one neglects, as we always do, a weak logarithmic dependence on 
p 2 induced by anomalous dimensions of the operators involved.) The third type of terms have 
been taken into account in eq. (7) in the leading order with a small total effect on (B — 1). 
Thus, all the contributions are under control within the approach and prove to be small. 

However, there exists another source of uncertainty - the use of the procedure of vacuum 
saturation to evaluate the v.e.v.'s of four-quark operators, which might (as simple estimates 
show) violate drastically eq. (9). Let us discuss the issue in some detail. 

To begin with, it is easy to show that the peculiar combination of four-quark operators 
entering eq. (7) transforms as a member of a 27-plet with respect to the (flavour) SU^(3) 
group. This means that the corresponding v.e.v's might only appear in the second order of the 
expansion in the strange-quark mass. It can be checked that this suppression is still operative 
if one takes into account next-to-leading corrections to the coefficient functions of the operators 
under consideration. However, it is interesting to get a direct estimate of the accuracy of the 
vacuum saturation hypothesis. To this end it is convenient to employ the QCD sum rules 
method in configuration space (x-space). It will be seen below that such x-space sum rules 
provide the unique possibility of direct estimation of the v.e.v.'s of four-quark operators. 

Indeed, let us consider the correlator 

{Tf(x)U0)) = n(-* 2 ), 

where j 1 * = u'y^d is the interpolating current of the p-meson. At x 2 — > the following x-space 
operator expansion holds: 



'2, 4 13 



6 ( L + 2 lE + I ) (u l5l x t a dd l5l x t a u)+ (10) 



2 



+ (-L + - lE - -) (^H^u^u + d^t a d)) 



+ o(l) 



where L = ln(— fi 2 x 2 /4), 7^ = 0.577. . ., and fi is the normalization point of the MS-scheme. 
Expansion (10) can be obtained from the corresponding one in p-space. Note that the val- 
ues of the constants added to L in the right-hand side of eq. (10) are fixed by the recipe for 
renormalization of the operator j^j^. Note also that expansion (10) remains valid after substi- 
tutions j M — > j'5 1 , 7 5 7a — > 7a, which correspond to considering the axial current interpolating 
the Ax-meson. 
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On the other hand the following Kallen-Lehmann representation takes place: 



U(-x 2 ) = p(s)A c (-x 2 ,s)ds, A c (x 2 ,s) = — —Vx^KUVxh;), (11) 

where the spectral density can be roughly approximated as follows: 

p(s) = F5(s — m 2 ) + a6(s — s ), 
the parameters F, m, sq, a being known from the analysis of the same correlator in p-space 

B 

The main advantage of relations (10) and (11) over similar sum rules in p-space comes 
from the fact that the contribution of the operator j^j^ on the right-hand side of eq. (10) 
is parametrically increased due to the contact term pictured in fig. 1 and corresponding to a 
disconnected diagram. This allows one to use eqs. (10) and (11) for determining the v.e.v. 
0'%) directly. 




x 



Fig. 1 - The origin of the contact terms in the x-space operator product expansion. 

An analysis of sum rules in x-space leads to the relations (see fig. 2 and 3; we took \x = 
1 GeV, a s (lGeV) = 0.3, (qq) = (-0.23 GeV) 3 ) 

(u^dd^u)/(qq) 2 = --(0.90 ± 0.15) , 
(uY^ddYl 5 u)/(qq) 2 = ^(0.84 ± 0.20) , 
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Fig. 2 - The curve showing the dependence of the numerical value of the v.e.v. 
(u'y ll dd'j ll u) on the variable \x\ — \/—x 2 . 




Fig. 3 - The curve showing the dependence of the numerical value of the v.e.v. 
(u^j^dd'y^j^u) on the variable \x\ = \J —x 1 . 



which are in agreement with the result of the vacuum saturation within an accuracy 30%. The 
relative discrepancy between v.e.v.'s of four-quark operators of VV and AA types is about 
10%; this fact indicates in favor of a small contribution of four-quark operators entering the 
right-hand side of eq. (7) to the value of B. 



3. - Let us discuss another calculation of B [10] employing the finite energy sum rules for 
the two-point correlator 

P(Q 2 ) = ijdx exp[iga;](0|Td(a;)6(0)|0) 

together with the method of effective chiral Lagrangians, the latter being used to fix the func- 
tional dependence of the corresponding spectral density on energy. The obtained value of B 
differs noticeably from both result (9) and the prediction of the vacuum saturation hypothesis. 
It should be noted that this calculation involved a large-size correction due to the operators 
m 2 , m 4 , mqq. This fact forced the authors to use a rather large "duality interval" Sq = 8 GeV 2 



in handling their sum rules. On the other hand, within the approach of paper [12] all these 



contributions are exactly summed up in the factorized term. Note also that the use of the 
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chiral perturbation theory assumes the independence of the form factor from the energy of 
the K° — K° pair, which can hardly be a reasonable approximation at energies as large as 



y/So = (2.5 -T- 3) GeV. The technicalities aside, we feel that the result of Pich and Rafael may 
well be somewhat underestimated. 

4. - To conclude, we have used the x-space sum rules to get a direct estimation of the 
accuracy of the vacuum saturation technique for operators ii'y^ddj^u and uj^j^ddj^^u. Our 
result is that the v.e.v.'s of these operators are equal to each other with an accuracy of about 
10%. This observation gives an additional support of the result of the calculation of the 
parameter B which was discussed in sect. 2. 
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